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co 
O 

o 

(N 
O 

Q 

(N 

(N 
> 
co 

(N 



CO 

O 

Ok 

cd 



X 



T. SKWARNICKI 

Department of Physics, Syracuse University, NY 13244, USA 
E-mail: tomasz@physics.syr.edu 

I review heavy quarkonium physics in view of recent experimental results. In particular, I discuss new results on 
spin singlet states, photon and hadronic transitions, D— states and discovery of the yet unexplained narrow X (3872) 
statě. 



1 Quarkonia 

Quarkonium is a bound statě of a quark and its anti- 
quark. Some properties of light and heavy quarkonia 
are compared to properties of positronium in Table 1. 
Unlike positronium, light quarkonia are highly rela- 
tivistic. They also contain mixtures of quarks of dif- 
ferent flavor and fall apart easily into other mesons. 
Charmonium (cč) was the first heavy quarkonium 
discovered and is less relativistic; the number of long- 
lived states below the dissociation energy (i.e. the 
threshold for decay to DĎ meson pairs) equals the 
number of long-lived positronium states. Bottomo- 
nium (66) is even more non-relativistic and has a 
larger number of long-lived states. The toponium 
systém would have been completely non-relativistic. 
However, weak decays of the top quark will domi- 
nate over the strong binding and long-lived states 
will not be formed. Therefore, charmonium and bot- 
tomonium play a speciál role in probing strong in- 
teractions. 

The states below open flavor threshold live long 
enough for electromagnetic transitions between vari- 
ous excitations to occur. The electromagnetic transi- 
tions compete with transitions mediated by the emis- 
sion of soft gluons. The latter materialize as light 
hadrons. Eventually the heavy quarks must anni- 
hilate into two or three hard gluons. Properties of 
these bound states and their decays are good testing 
grounds for QCD in both the non-perturbative and 
perturbative regimes. 

The first heavy quarkonium bound statě above 
the DĎ or BB threshold acts as a factory of heavy- 
light mesons. The heavy quarks trapped in these 
mesons ultimately decay via weak interactions. This 
is a good pláce to look for physics beyond the stan- 
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dard model as discussed by Y. Grossman at this 
Symposium. 1 Many measurements of electroweak pa- 
rameters are obscured by strong interactions. This 
provides an important motivation for trying to un- 
derstand details of strong interaction phcnomcna. 
P. Lepage pointed out that at least some of the 
new interactions to be discovered are likely to be 
strongly coupled, further motivating detailed studies 
of QCD. 2 

Heavy quarkonia offer two small parameters: ve- 
locities (v) of constituent quarks and the strong cou- 
pling constant (a s ) in annihilation and production 
processes. The expansion of the full theory in these 
parameters allows for effective theories of strong in- 
teractions: in the past - purely phenomenological 
potential models; more recently - NRQCD. 3 Lattice 
QCD calculations are also easier for heavy quarkonia 
than for light hadrons. 

2 Hadroproduction 

Annihilation of n 3 Si~- {4>, T) states to /i + ^t _ and to 
a lesser extent e + e~ makes it possible for experimen- 
talists to fish out these states from large backgrounds 
in hadroproduction experiments. This is how these 
states were (co)discovered! There is also some access 
to the production of the excited states by addition 
of a transition photon or 7r + 7r~ pair. So far, except 
for the initial discovery, the hadroproduction exper- 
iments have not played an important role in spec- 
troscopy or decay studies, but have made interest- 
ing production measurements. Heavy quarkonium is 
also a useful probe for determining the structure of 
the target (e.g. to probe for gluon content or the 
presence of quark-gluon plasma) . 

Non-relativistic QCD (NRQCD) is a favorite 
theoretical approach used to describe the produc- 
tion data. 4 Various diagrams are ordered according 
to powers of a s and v. The diagrams can also be 
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classified into color-singlet processes, in which the 
heavy QQ pair is produced colorless and therefore 
can directly form a bound statě, and color-octet pro- 
cesses, in which the QQ pair has color which must 
be emitted away by soft gluon (i.c. long distance) in- 
teractions. Initial attempts to deseribe charmonium 
production at the Tevatron with only color-singlet 
processes (CS model) 5 failed spectacularly, revealing 
the importance of color-octet diagrams. However, 
the color-evaporation model (CEM), 6 which allows 
color-octet processes but neglects their ranking in v, 
also fits the data well. All theoretical approaches in- 
volve free parameters, which are fixed by fits to the 
data, diminishing their effectivc diffcrcnccs. 

Polarization of produced heavy quarkonia should 
lead to a better diserimination between NRQCD and 
CEM. NRQCD predicts an inerease in polarization 
with higher p tl while CEM predicts no polarization 
at all. The present data do not provide any evi- 
dence for inerease in polarization with p t , however 
the experimental errors are too large to draw any 
firm conclusions. 

Thcre is a large range of kinematical regimes 
and differential cross sections for charmonium pro- 
duction studies at HERA. Both Hl and ZEUS con- 
tributed a number of papers on this topič to this 
Symposium. 7 NRQCD predicts that matrix elements 
should bc univcrsal. However, it is dimcult to recon- 
cile all Tevatron and HERA data with a consistent 
set of NRQCD matrix elements. It is quite possible 
that the charm quark is just not heavy enough for 
the NRQCD approximation to work well. Unfortu- 
nately, data on bottomonium production are either 
non-existent or have large experimental errors. 

For a more complete review of this topič see e.g. 



Krámer. 

3 Clean Production Environments 

Most of what we know about heavy quarkonium 
states and their decays comes from experiments at 
clean production environments, which are time re- 
versals of simple decay modes (see Fig. 1). 

Vector states (J PC = 1 : n 3 Si, n 3 Di) decay 
to lepton pairs and thus can be directly formed in 
e + e~ collisions. Production rates are large compared 
to the other e + e~ processes, thus the backgrounds 
are small and these states can be studied in both in- 
clusive and exclusive decay modes. Dedicated runs 
are needed for cach vector statc. The other states 
can be reached via photon and hadronic transitions, 
however, their scope is limited by the transition sc- 
lection rules and branching ratios. 

Spin and 2 states with positive C-parity 
(J PC = (T+,0+ + ,2++: r^So, n 3 P {)a ) decay to two- 
photons, thus can be formed in two-photon collisions 
at the e + e~ colliders. No dedicated runs are needed 
sincc the two-photon flux populates a wide range of 
possible quarkonia masses. However, the fiux drops 
quickly with the energy and production rates for 
heavy quarkonia are small compared to the dominant 
e + e~ processes. To combat these large backgrounds 
the quarkonia states must be detected in simple ex- 
clusive decay modes. So far this formation method 
has succeeded only for charmonium states. 

States of any quantum numbers can decay to pp 
via two- or three-gluon annihilation, however, such 
couplings are small. Low energy p beams annihilat- 
ing with proton-jet targets were successfully used to 
form charmonium states. The beam energy must be 
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tuned to form a specific quarkonium statě. Since the 
background cross-sections are very large, the quarko- 
nia must be detected in simple exclusive finál states 
to take advantage of the highly constrained kinemat- 
ics. 

Decays of B mesons also offer a clean production 
environment for charmonium states. Again states of 
any quantum number can be formed. Since the pro- 
duction rates are rather small, experimentalists must 
restrict themselves to specific exclusive finál states to 
take advantage of the B mass constraint (and beam 
energy constraint if produced at the T(45) resonance 
in e + e~ collisions). 



e+ e + 




q q 

Figuře 1. Clean production environments for heavy quarkonia 
together with corresponding decay modes for most of them. 

4 Results from pp Experiments 

A H2 gas jet is used as a target for an antiproton 
beam in the p accumulator ring. This technique was 
pioneered at the CERN-ISR by the R704 experiment 
(1983-84). 8 It was later used at Fermilab with the 
E760 detector (1989-91) 9 followed by its major up- 
grade E835 (1996-97, 2000). Since this was a non- 
magnetic detector, it was limited to specific finál 
states containing electrons and photons. Charmo- 
nium states are observed directly as resonance peaks 
in the signál event yield measured as function of the 
center-of-mass energy. Masses and widths of the ob- 
served states can be measured with high precision, 
since the accuracy depends on knowledge of the beam 
energy rather than resolution of the detector. The 
Fermilab experiments determined masses and widths 



of the Xcj states with unprecedented accuracy. 9 ' 10 
They also studied radiative transitions from these 
states. 11 The most recent results include detection of 
the XcO resonance in decays to interfering with 

continuum production of this finál statě 12 and obser- 
vation of the singlet ^(l 1 ^) statě in annihilation to 
two photons. 13 The radiál excitation of the latter was 
not found in scans by E760 and E835. 14 Only the re- 
gion around the mass reported by the Crystal Ball 
experiment 15 was scanned. 

While the singlet ?7 c (l 1 S'o) was observed in many 
different environments, the only hints for the sin- 
glet /i c (l 1 Pi) statě come from the pp experiments. 
R704 provided some inconclusive evidence for this 
statě in 1984. 16 Better evidence for this statě, with 
the mass close to the center-of-gravity mass of the 
X c (l 3 Pj) states, was reported by E760 17 via the de- 
cay chain: h c — > 7r°J/-0, J/ip — ► e + e~. E835 took 
more scans in this mass range, with larger statistics 
and an improved detector. Recently, rumors of ab- 
sence of the h c statě in the E835 data was reported 
in print by several non-E835 authors 18 However, the 
official word from the E835 collaboration 19 is that 
they are analyzing all available decay channels and 
are not ready to make any definite statements yet. 

Average: (2979.9 ±1.0) MeV 
CL=0.5% Scale Factor=1.5 
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Figuře 2. Measurements of the r) c mass. The new measure- 
ments (BES 03, 21 BABAR 03, 22 E835 03, 13 BELLE 02 23 and 
CLEO-c 03 24 ) are indicated by a star. See PDG 2002 20 for 
the rcfcrcnccs to the older measurements displayed herc. The 
vertical bars indicate the average value (solid) and its error 
(dashed). The error includes the scale factor. The thick hor- 
izontál bars on the right give the relative weight of each ex- 
periment into the average value. 
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5 Charmonium Singlet States 

There are many new measurements related to the 
charmonium singlet states, in addition to the results 
from the pp experiments mentioned in the previous 
section. For example, there are 5 new measurements 
of the rjci^So) mass 13 ' 21 " 24 (see Fig. 2). BES-II, 21 
which reconstructed the r\ c in 5 different decay modcs 
in the J/ip data, achieved the smallest overall er- 
ror. Consistency of the old and new measurements 
is questionable. A new prcliminary measurement of 
the r\ c width by BABAR 22 constitutes an even big- 
ger experimental puzzle (see Figs. 3-4). The mea- 
surement error claimed by BABAR is much smaller 
than in any other measurement. At the same time 
the centrál value, (33.3 ±2.5 ±0.5) MeV, is the high- 
est ever measured; completely inconsistent with the 
previous world average value 20 (16.0Í 3 ; 2 ) MeV. 
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Figuře 3. Observation of 77 — > r] c (í8) and 77 — > ri c (2S) by 
the BABAR experiment. 22 A peak due to e + e — — > ~(J/ip is 
also visible. These data serve as a precision determination of 
the r/ c (lS) width and confirmation of the r] c (2S) statě at its 
new mass value. 

There have also been dramatic developments 
concerning the 'q c (2 1 So) statě. The Crystal Ball ex- 
periment at SPEAR claimed observation of this statě 
over 20 years ago. 15 They observed a peak in the 
inclusivc photon energy spectrum from ip(2 3 Si) de- 
cays, which they interpreted as a magnetic dipóle 
photon transition, ip(2S) -> ~/r] c (2S). 15 The Crystal 
Ball result corresponded to the hyperfine mass split- 
ting in the radiál excitation, which was only slightly 
smaller than in the ground statě (92 MeV vs. 117 
MeV). Last year BELLE observed both r] c states in 
B -> Kr) c (nS), r) c {nS) -> K s K+tt-. 25 The r] c (2S) 
statě appeared at much higher mass than in the 
Crystal Ball measurement, thus reducing the corre- 
sponding hyperfine splitting. BELLE also observed 
both r] c states in the spectrum of the mass recoil- 
ing against J/ip in continuum e + e _ annihilation to 
J/i/jX? 6 These results have been updated this year 
with larger statistics. 27 The r] c (2S) mass obtained 



from these data differs by 1.9 standard deviations 
from the other BELLE result, but is stih significantly 
higher than the Crystal Ball value. 

Average: (25.0 ± 3.3) MeV 
CL=0.05% Scale Factor=1.8 
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Figuře 4. Measurements of the T) c width. See Fig. 2 for refer- 
ences and explanation. 

The production rate for charmonium pairs ob- 
served by BELLE in the latter analysis is surpris- 
ingly high. They also observe a large e + e~ anni- 
hilation rate to J/tp plus open charm. They claim 
(82 ±15 ±14)% of prompt J/tp production [Pj/^ > 2 
GeV) at 10.6 GeV center-of-mass energy is associ- 
ated with another cč pair in the event. 27 This rate 
is an order of magnitude higher than that predicted 
theoretically, and constitutes the biggest puzzle in 
quarkonium production physics. 28 

This year the BABAR and CLEO experiments 
confirmed the rj c (2S) at its heavier mass by for- 
mation in the two-photon collision process 22 ' 29 (see 
Fig. 3) . The discrepancy with the Crystal Ball mass 
measurement is resolved by the CLEO-c experiment, 
which has been able to remeasure the inclusive pho- 
ton spectrum from the ip(2S) for the first time since 
the Crystal Ball experiment (see the next section). 
While CLEO-c agrees well with the Crystal Ball on 
all other El and Ml transitions from the ip(2S), the 
direct Ml transition to the rj c (2S) is not observed 
at the photon energy claimed by the Crystal Ball. 
The upper limit on the rate for this transition set 
by CLEO-c disagrees with the rate measured by the 
Crystal Ball. Therefore, the peak observed by Crys- 
tal Ball could not be due to rj c (2S) production. 

Measurements of the r] c (2S) mass are summa- 
rized graphically in Fig. 5. All new mass mea- 
surements are consistent with each other. The to- 
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tal width of this statě, based on the BELLE 25 and 
BAB AR 22 results, is (19 ± 10) MeV. 

The new mass of the rj c (2S) statě decreases the 
hyperfine mass splitting for the 2S'-states by a factor 
of 2 compared to the old value. Since a wrong value 
prevailed for 20 years, it is interesting to check for an 
experimental bias on the phenomenological predic- 
tions for this splitting. A sample of potential model 
predictions for the ratio of the hyperfine mass split- 
ting, R HF = (M 4 , { 2s ) ~M ric (2S))/( M ^(is)- M Vc (iS)), 
is compared to the old (R = 0.79) and the new 
(Rhf = 0.412±0.028) experimental values in Fig. 6. 
It appears that many old calculations were stretched 
to accommodate the old result, though there were 
a few that had the courage to contradict the data. 
The ratio of hyperfine mass splitting can be related 
to the ratio of the leptonic widths of the triplet states 
(r ee ) using perturbative QCD. The new value agrees 
well with the prediction by Badalian and Bakker, 31 
Rhf = (0.48 ± 0.07), based on this approach. In- 
stead of using the measured values of T ee , Reck- 
siegel and Sumino 32 extended the use of perturbative 
QCD to the extraction of the interquark potential 
at short distances relevant for the spin-spin forces. 
Their prediction for the hyperfine mass splitting ra- 
tio, Rhf = 0.42, is in good agreement with the data, 
but the absolutc values for the mass splittings are un- 
derestimated. Lattice QCD calculations 33 predicted 
Rhf — 0.5, not far from the measured value. 

Average: (3637.7 ±4.4) MeV 
CL=14% Scale Factor=1.3 



3620 3640 
Mass (MeV) 



Figuře 5. Measurements of the r? c (25) mass (see the text for 
the references). The thick horizontál bars on the right give the 
relative weight of each experiment into the average value. The 
Crystal Ball measurement was excluded from the average. 
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Figuře 6. Measurements and predictions for the ratio of the 
hyperfine mass splittings in the 25 and 15 states in charmo- 
nium. The old and new experimental values are indicated by 
vertical bars. The dashed lineš around the new value show 
the experimental error. The histogram shows predictions of 
a sample of potential models. 30 Only models which gave a 
good fit to the masses of all known charmonium states were 
included. The shaded part of the histogram represents the 
models published before the new measurements of the t; c (25) 
mass. Values predicted by scaling from the measured leptonic 
widths ratio 31 (r ee ), by perturbative QCD alone 32 (pQCD) 
and by lattice QCD 33 (laQCD) are also indicated. 

6 Radiative Transitions 

Since sizes of heavy quarkonia are small or compa- 
rable to the wavelengths of transition photons, se- 
lection rules are given by multipole expansion. As 
the systém is non-relativistic, electric dipóle tran- 
sitions (AL = 1, AS* = 0) are much stronger than 
magnetic dipóle transitions (AL — 0, AS* = 1). 
The latest improvements in measurements of these 
transitions come mostly from the CLEO experi- 
ment, which is equipped with an excellent CsI(Tl) 
calorimeter. Recently, after increasing statistics of 
the T(l,2,35) data samplcs by an order of mag- 
nitude, CLEO turned the beam energy down and 
collected also ip(2S) data. This was the beginning 
of the CLEO-c phase. Inclusive photon spectra col- 
lected at iji(2S) and T(2S) by CLEO are compared 
in Fig. 7. The dominant peaks are due to El tran- 
sitions, 2 3 Si — » 7l 3 Pj. The peaks in bottomo- 
nium systém appear smaller because of the increased 
7r — > 77 background induced by higher particle mul- 
tiplicities. The peaks are more crowded together re- 
flecting the smaller fine structure of the l 3 Pj states 
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in the more non-relativistic T systém. The El cas- 
cade lineš, l 3 Pj — ► jl 3 Si are also visible. The char- 
monium data also reveal a small peak due to the 
hindered Ml transition, 2 3 S'i -> 7l 1 S , . This is the 
first confirmation of this transition since the origi- 
nál observation by the Crystal Ball. As discussed 
in the previous section, the direct Ml transition, 
2 3 S'i — > 72 1 S'o, is ruled out for the photon energy, 
r] c (2S) width and rate claimed by the Crystal Ball. 
The preliminary CLEO-c results 45 for the ip(2S) pho- 
ton lineš are: 2?(V>(2 3 Si) -> IXÁ^Pj)) = {9.75 ± 
0.14± 1.17, 9.64±0.11±0.69,9.83±0.13±0.87}% for 
J = {2,1,0}, B(V(2 3 5i) -> ^vÁ^So)) = (0.278 ± 
0.033 ± 0.049)% and B{ip{2 3 S í ) -► 7?7 c (2 1 S'o)) < 
0.2% for £ 7 = (91 ± 5) MeV and r(í? c (25)) = 8 
MeV (90% CL). 
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Figuře 7. Inclusive photon spectrum from 2 3 Si decays in the 
cč (top) and bb (bottom) systems measured with the CLEO 
detector. The data correspond to about 1.5M i/>(2S) and 9M 
T(2S) decays. 

None of Ml transitions, hindered or direct, are 
observed for the bottomonium. This is not a big 
surprise, since the backgrounds are much higher and 
the expected branching ratios are smaller in the more 
non-relativistic bb systém. Searches for the singlet rfo 
statě in two-photon collisions at LEP has also yielded 



upper limits only, thus no singlet states have been 
observed in bottomonium so far. The Ml rate mea- 
surements can be translated into values of the cor- 
responding matrix elements, which are compared to 
theoretical predictions in Fig. 8. The matrix ele- 
ments for the hindered Ml transitions are expected 
to be very small, since they are generated by rela- 
tivistic and finite size corrections. Thcrefore, they 
are difficult to predict. Only very recent calcula- 
tions are consistent with all charmonium and bot- 
tomonium data. Similar comparisons for the El ma- 
trix elements is shown in Fig. 9. Non-relativistic cal- 
culations (hollow circles) overestimate the El rates 
in charmonium. The predictions with relativistic 
corrections (filled triangles) can reproduce the data. 
Relativistic effects in the dominant El transitions 
in bottomonium are much smaller. However, rela- 
tivistic calculations are needed to reproduce the rare 
3 3 Si — > 7l 3 P j transition rate. This matrix element 
is small because of the cancellations in the integrál 
for the El operátor between the initial and the finál 
statě wave functions 39 for which the number of nodes 
differ by two. 
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Figuře 8. Measured and predicted values of matrix elements 
for Ml transitions in heavy quarkonia. The measured values 
are calculated from the world average branching ratios and to- 
tal widths of the triplet S states. The centrál value and error 
bars for the cč data are indicated with solid and dashed lineš 
respectively. For the bb data, allowed ranges from the prelim- 
inary CLEO analysis 45 are shaded (90% CL). The theoretical 
predictions 35 (points) are ordered according to the publication 
date. 
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Figuře 9. Measured and prcdictcd values of matrix clcmcnts 
for El transitions in heavy quarkonia averaged over different 
spins of the triplet P states. The measured values are calcu- 
lated from the world average radiative branching ratios and 
total widths of the triplet S states. The centrál values and 
error bars for the measured values are indicated with solid 
and dashed lineš respectively. Circles (triangles) show non- 
relativistic (relativistic) calculations. The relativistic calcu- 
lations are averaged over spins with the same weights as the 
data. The predictions 36 are ordered according to the publica- 
tion date. 

7 Hadronic Transitions 

Heavy quarkonia can also change excitation levels by 
emission of soft gluons, that turn into light hadrons. 
The multipole expansion approach has proven to be 
useful also for hadronic transitions 37,38 explaining 
their gross characteristics. The most prominent are 
7T7r transitions among n 3 Si states, which can be me- 
diated by two-gluon emission of the El -El type. In 
fact, these are dominant decays for both tp(2S) and 
T(2S). The ratio of the measured widths for the 
these transitions agrees with a suppression by about 
a factor of 10 predicted by the multipole expansion 
model 38 due to the smaller size of the bb systém. The 
multipole expansion model is also able to explain the 
M (tttt) distributions for these transitions, which fol- 
low the same pattern and peak at high values. Dip- 
ion transitions from T(3S) to T(2S) and T(15) are 
observed too. Thcy have smaller rates either due to 
the phase space suppressions or cancellations in the 
dipóle matrix element integrál for the 3S — ► mrlS 
transition. The M(irir) distribution peaks at low and 
high mass values for the latter transition, which re- 
veals some dynamics beyond the multipole expansion 
approach. In charmonium, an rj transition has been 



observed between the triplet 2S and 15 states. 20 This 
transition is of a magnetic type (E1-M2), thus it has 
a smaller rate than the irir transition. A ir° tran- 
sition is also observed, but with a tiny rate due to 
isospin violation. Nonc of these transitions are ob- 
served between the T states, 45 which is not surprising 
since there is an additional suppression by at least a 
factor 1/rriQ for this type of transition. 
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Figuře 10. Recoil mass against the lowest energy photons in 
four-photon cascade between the T(3S) and T(IS), consis- 
tent with X(,(2P2 J i) and Xb(lP2,i) as intermediate states, as 
observed by CLEO. The cascades via the T(2S) statě are sup- 
pressed, thus the remaining events are T(1D) candidates. The 
recoil mass is a measure of the T(1D) mass. The fit to the 
data (solid line) assumes production of just one T(1D) statě, 
with a natural width much smaller than the experimental res- 
olution. The observed events are most likely dominated by 
the production of the T(1D2) statě. The preliminary CLEO 
results give (10161.1 ± 0.6 ± 1.6) MeV for the mass of this 
statě. 

The transitions among triplet nS states dis- 
cussed above were hrst observed over 20 years ago. 
A large number of other hadronic transitions are 
possible, 39 especially for the bb systém which has a 
large number of long-lived states. Such transitions 
had not been observed until recently, when CLEO re- 
ported the hrst observation of Xb(2 3 P 2 ,i) — > wT(lS) 
transitions. 40 The phase-space for these transitions 
is so small, that this decay is forbidden for the 
Xb(2 3 Po). The measured branching ratios for the 
Xb{2 3 Pi) and x&(2 3 Í2) states are of the order of a 
couple of percent ((1.6 ± 0.3 ± 0.2)% and (1.1 ± 0.3 ± 
0.1)% respectively), in špite of the phase space sup- 
pression, which reveals that the underlying transition 
is quitc strong. In fact, this transition is of chromo- 
electric type, El -El -El (three gluons are needed to 
generate it). Voloshin 41 pointed out that since the 
matrix element does not depend on the spin of the 
2Pj statě, transition branching ratios for J = 1 and 
J = 2 should be comparable (as the phase space 
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factors approximately cancel the effect duc to the 
smaller total width of the J = 1 statě). The data 
are consistent with this prediction. 

As hadronization probabilities are difřicult to es- 
timate, uncertainties in absolute rate predictions for 
hadronic transitions are usually very large. When 
rcsonances dominate the transition, there are often 
no theoretical predictions for the rate. For example, 
there are no rate estimates for the u> transition dis- 
cussed above. Predictions for other yet unobserved 
types of transitions vary by orders of magnitude. The 
predictions based on a model developed by Yan 38 
and Kuang 42 tend to be much larger than the pre- 
dictions based on a differcnt approach introduced by 
Voloshin, Zakharov, Novikov and Shifman. 43 
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o(AM C0fl ) (MeV) 

Figuře 11. The difference between the predicted and the mea- 
sured T(1 3 D2) mass (vertical axis) versus overall quality of 
the potential model expressed as the r.m.s. of the differences 
for the center-of-gravity masses of long-lived bb states (hori- 
zontál axis). Predictions for the 1S, 2S, 3S, IP, 2P and ID 
masses are included in the r.m.s. calculation for all models 
displayed here. 46 
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Figuře 12. Scarch for Y(l 3 r>2) — > 7r+7r _ T(lS') transitions. 
Recoil mass against two photons in ■yjir^n^ l + l~ events is 
plotted. The expected signál position is indicated by the ar- 
row. 

One suitable pláce to test various predictions are 
dipion transitions from l 3 Z?j to l 3 >Si. Last year 



the T(l 3 Dj) states were discovered by CLEO in 
the four-photon cascade: 44 T(3 3 ,Si) -> 7Xb(2 3 P/ 2 ), 
X6 (2 3 P /2 ) - 7 T(1 3 A/ D ), T(1 3 A/J - 7X6(1^), 
Xb{l 3 P. h ) -> 7T(1 3 5!) followed by T^S^ -» 
CLEO has analyzed more data since then. The re- 
coil mass distribution against the hrst two photons 
in the cascade is consistent with the detector reso- 
lution (Fig. 10). Thus the data are dominated by 
production of just one T(í 3 Dj D ) statě. Theoreti- 
cally, production of the Jd = 2 statě is expected to 
dominate. This spin assignmcnt is also favored by 
the experimental data. 44 The measured mass fits the 
predictions of the potential models well, especially 
those which provide a good fit to the other narrow 
bb states. This is illustrated in Fig. 11. The lattice 
QCD calculations also reproduce the observed mass 
very well. 2 CLEO determines 45 the four-photon cas- 
cade branching ratio to bc (2.6 ± 0.5 ± 0.5) x 10~ 5 . 
Within the errors this rate is consistent with the 
predictions by Godfrcy and Rosner 47 3.8 x 10~ 5 
(2.6 x 10~ 5 for Jd = 2 alone), based on the El ma- 
trix elements and estimates of hadronic widths of the 
T(1 3 L>,/J and Xfc(2 3 P/ 2 ) states. CLEO also looked 
for T(1 3 A/J Tr+Tr-r^S*!) replacing the third 
and fourth photons in the cascade with a 7r + 7r~ pair. 
The recoil mass against the remaining two photons 
in the cascade is plotted in Fig. 12. Since no signál 
is found, upper limits on the product of branching 
ratios are set. They are compared in Table 2 to the 
theoretical predictions calculated by Rosner 48 using 
predictions for ?(1 3 D Jd ) tt+tt-T (l 3 S*i) width 
by various authors 42 ' 49 ' 50 together with the El- 
photon matrix elements and estimates of hadronic 
widths of the T(l 3 Z)j D ) and Xb(2 3 Pj 2 ) states. Large 
T(l 3 Dj D ) 7r+7r-T(l 3 Si) widths, as predicted 
by Kuang- Yan, 42 are ruled out. Better experimen- 
tal sensitivity is needed to test the other models. 
Voloshin pointed out 51 that T(l 3 DjJ ??T(l 3 5i) 
transition may be of comparable strength to the irw 
transition between these states. CLEO doesn't find 
this transition either and sets a prcliminary upper 
limit of ^(T^!) - 1X b(^P.h)) B( Xb (^Pj 2 ) - 
7 T(1 3 A/J) B(T(l 3 Dj D ) - ryT(l 3 Si)) < 2.3xl0- 4 
at 90 % CL. 

There are also new results on l 3 Di — > n + Tr~l 3 Si 
transitions in the charmonium systém. Here, the- 
oretical situation is complicated by mixing of the 
?/>(l 3 -Di) statě with the V'(2 3 5'i). The observed statě, 
?A(3770), is also above the open flavor threshold, 



Table 2. Theoretical predictions and prcliminary CLEO up- 
per limits (at 90% CL) on 0(T(3 3 Si) -> 7Xf,(2 3 Pj 2 )) 8(x b 
(2 3 P /2 ) ^jT(1*Dj d )) B(TQ.»Dj D ) -» 7r+7r-T(l35i)) in 
units of 10 4 . The theoretical predictions come from the pá- 
per by Rosner. 48 The first row corresponds to the cascadc 
via the Jjj = 2 statě observed in the four-photon cascade by 
CLEO. The second row corresponds to any T(l 3 Dj r) ) statě 
with mass in the 10140 - 10180 MeV interval. 
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therefore it has a large width for decay to DĎ which 
suppresses branching ratios for any transitions to 
the other charmonium states. BES claims observa- 
tion of such a transition in the data consisting of 
5.7 x 10 4 ^(3770) decays. 52 Their reconstruction ef- 
ficiency is 17%. They observe 9 signál candidates 
with an estimated background of 2.2 ± 0.4 cvents 
(see Fig. 13a). The corresponding branching ratio is 
(0.59 ± 0.26 ± 0.16)%. Such a large branching ratio 
would favor the Kuang-Yan model 53 contrary to the 
66 results discussed above. CLEO-c has also analyzed 
their first V (3770) sample. 45 They have a smaller 
sample of ^(3770) decays (4.5 x 10 4 ) but larger recon- 
struction efficiency (37%). They have only 2 events 
in the signál region, consistent with the estimated 
background level (see Fig. 13b). They set a prclim- 
inary upper limit S(^(3770) -> tt+tt- J/ip) < 0.26% 
(90% CL). This result is not inconsistent with the 
BES value because of the large experimental errors 
in both measurements. However, the CLEO result 
indicates that it is premature to favor the Kuang- 
Yan model on the basis of the V- (3770) data. BES 
is already analyzing a larger data sample and CLEO 
is expected to inerease their statistics by an order of 
magnitude this fall, thus more aceurate results are 
expected next year. 

8 New Particle Discovered by BELLE 

Decays to tt+tt- Jj ip are also a way to search for other 
charmonium states, for cxamplc ft, c (l 1 Pi), ip(l 3 D 2 ), 
ip(l 1 D2) etc. Many of these states cannot be di- 
rectly formed in e+e~ annihilation, however they can 
be produced in pp annihilation, B decays or fragmen- 
tation processes. BELLE inspected resonance struc- 
tures in the tt + tt~ J /ip systém produced in the decay 
B ± -+ (tt+tt" J/VO, J/ip -» l+r. The distri- 
bution of AM = M (tt+tt' J/ip) - M(J/ip) observed 
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Figuře 13. Searchcs for «/)(3770) — > tt+it~ J/tp transitions. 
The expected signál positions are indicated with arrows. The 
dominant peak in cach distribution is duc to e+e~ — > ytp(2S), 
ip(2S) -> w+ir-J/ip. (a) BES results. M(l+l~) mass for kine- 
matically constrained 7r+7r~Z+Z — events is plotted. The signál 
is expected at the J/ip mass. (b) CLEO results. The difference 
between the center-of-mass energy and the recoil mass against 
the tt + it~ pair is plotted with no kinematic constraints, but 
after the eut on M(l + l~) around the J/^ mass. The signál is 
expected at the mass difference between the t/>(3770) and the 
<//</>■ 

by BELLE 54 in a sample of 3 x 10 8 B mesons is 
shown in Fig. 14. The prominent peak is due to the 
tp(2 3 Si). There is also a smaller but very significant 
peak observed at a higher mass. By performing an 
unbinned maximum likelihood fit to AM, the beam- 
constraint B meson mass and energy difference be- 
tween the B candidate and the beam energy BELLE 
finds 35.7 ± 6.8 events in the second peak with a 
statistical significance of 10.3 standard deviations. 
The mass determined from the AM peak position 
is (3872.0 ± 0.6 ± 0.5) MeV. The observed width of 
the peak is consistent with the detector resolution, 
therefore the new statě is long-lived (r tot < 2.3 MeV 
at 90% CL). 

The invariant mass of the 7r + 7r~ systém for the 
signál events is strongly peaked at high mass values. 
The peaking is stronger than predicted by Yan 38 with 
the multipole expansion approach for S — > S transi- 
tions, and much stronger than predicted for D — > S 
transitions in the quarkonium systém, as illustrated 
in Fig. 15. The observed dipion mass distribution 
is suggestive of the isospin violating Y(3872) — > 
p°J/ip process instead. Isospin violation can be ex- 
perimentally verified by measuring T(Y(3872) — > 
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7r°7r°J/V)/r(X(3872) -> tt+tt- J/ip). 55 ' 56 For isospin 
conserving 7r7r transitions, this ratio should be ap- 
proximately 1/2, whereas p° does not decay to 7r°7r°. 

Soon after BELLE's announcement at the 
Lepton-Photon Symposium, the V (3872) particle 
was confirmed by CDF in 7T + 7r~ J/tjj, J/tp — > ^í + /í~. 57 
In 220 pb" 1 of Run-II data they observe 704 ± 67 
signál events. The preliminary mass measurement, 
3871.4 ± 0.7 ± 0.4 MeV, is consistent with BELLE's 
result. Their data also show peaking at high dipion 
mass. 
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Figuře 14. Distribution of M (tt+ tt~ 1+ 1~ ) - M(l+l~) as ob- 
served by BELLE for B decay candidates. 54 The first peak 
corresponds to the ip(2S) resonances. The second peak repre- 
scnts the X(3872) particle. 

Perhaps the most puzzling property of the newly 
discovered particle is its mass. Within the experi- 
mental errors it coincides with the sum of D° and 
D*° masses: m x - m D o - m D ,a = 0.2 ± 0.9 MeV 
(using the average of the BELLE and CDF rcsults). 
One cxciting intcrprctation, which actually predicts 
this to be the case, is a DĎ* (+ĎD*) molecule. 
Constituent mesons would likely be in a relative S- 
wave creating a J p = 1 + statě. "Molecular char- 
monium" was first discussed in the literatuře in the 
mid-seventies 58 and was initially introduced to ex- 
plain the complicated coupling of the e + e~ reso- 
nances above the open flavor threshold to various 
decay modes involving D and D* meson pairs. A 
satisfactory description of the e + e~ data was later 
achieved within a simple cč bound statě model. 59 ' 60 
Howevcr, since meson molecules were also proposed 
in the context of light hadron spectroscopy, 61 the 
concept of molecular charmonium did not go away. 
Interactions in the DĎ* systém (also BB*) were 
found to be attractive when described by the pion- 



exchange force. 62,63 No such force would exist in the 
DĎ (or BB) systém. Quantitative estimates showed 
that the DĎ* systém could be only loosely bound if 
at all, with better prospects for the BB* molecule. 
Tornqvist showed that in the limit of isospin sym- 
metry only an isoscalar molecule would be bound. 62 
Isospin is expected to be broken since the binding 
energy is comparable to the isospin mass splittings 
(the D+D-* threshold is 8 MeV above the D°Ď°* 
threshold) . Close and Page 64 argued that the isospin 
breaking does not change the conclusion that there 
is only one molecular systém expected near the DD* 
threshold. The loose binding makes it difficult for 
the molecule to rearrange the quark content in the 
meson subcomponents, thus strong decays to a char- 
monium statě plus light hadrons are expected to have 
small widths. Decays to -k + tt~ J/ip could proceed via 
the isospin violating p°J/ip channel consistent with 
the BELLE's dipion mass distribution. The molec- 
ular interpretation predicts that decays to D Ď°ir° 
and D°Ď°"f should oceur, since the D°* is almost 
on shcll. The ratio of the widths for these decays 
should be approximately 3:2 and their sum 60-100 
keV. 62 ' 64,65 ' 56 Therefore, the molecular model is con- 
sistent with the narrow width of the observed statě. 
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Figuře 15. Dipion mass distribution in BELLE's data 66 for 
X(3872) — * 7r+7r~ J/ip. The data are compared to the the- 
oretical shapes predicted by Yan 38 for D — > S (solid line) 
and S — » S (dashed line) transitions in the quarkonium sys- 
tém. The distributions were normalized to the same number 
of events. The efneiency dependence on M(tt + tt~) has been 
neglected in the comparison. 

Traditional charmonium states can also be nar- 
row at 3872 MeV if they cannot decay to DĎ. 
Possible candidates, ordered according to inereas- 
ing mass predictions are tp(l 3 D 2 --), r] C 2(l 1 D 2 -+), 
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Xc (2 3 P 1++ ), h c {2 1 P 1+ -) and ? 7c (3 1 5 - + ). 

The lower excitations of the r\ c are already 
broader than the upper limit on the natural width 
of the new statě, thus we can discard the r] c (3 1 S a -+) 
possibility. The other listed states can satisfy this 
upper limit. For the remaining options, radiative El 
transitions to the other charmonium states should 
dominate over irirJ/ip. This is especially true for 
the singlet states, for which the dipion transition 
would involve chromomagnetic interactions (spin- 
flip). Since the ft. c (l 1 P 1 +-) has not been seen in B 
decays, the /i c (2 1 P 1 +-) hypothesis is unlikely. 

BELLE has also just observed B+ — > 
#±^(3770) 0(3770) DĎ) for the first time. 67 
The measured branching ratio, (4.8±1.1±1.2) x 10 -4 , 
is comparablc to 20 B(B ± -» K+i){2S)) = (6.8 ± 
0.4) x 10~ 4 . Theorctically, B{B ± -> K ± ý(l 3 D 2 )) is 
expected to be 1.6 times larger than the production 
of V(3770) (assumed to be ^(l 3 L>i)). 68 Together, 
with BELLE's rcsult, 54 B(B ± -» K ± X)B(X -» 
tt+tt- J/VO = (0.063 ± 0.012 ± 0.007) xB(B ± -> 
K ± ^(2S))B(^(2S) -► tt+tt-J/V'), the ^(1 3 £> 2 ) in- 
terpretation requires B(ip(l 3 D 2 ) — > tt+tt- J/ip) to 
be about 1 to 3%, much smallcr than the naivě 
expectations 69 ' 56 for the dominant radiative decay of 
this statě, S(V>(1 3 -D 2 ) -> 7Xc(l 3 -Pi)) > 50%. Also, 
as discussed in the previous section the T(1 3 D 2 ) 
statě was observed via its radiative decay, whereas 
only an upper limit on the decay to 7r + 7r _ T ex- 
ists. However, BELLE observes no evidence for 
the photon transitions to Xc(l 3 -Pi) and sets the fol- 
lowing 90% CL limit: £(X(3872) -► 7X c (l 3 Pi))/ 
B(X(3872) -► tt+tt- J/^(l 3 5i)) < 0.89. Coupled 
channcl effects can be big (proximity of the DĎ* 
threshold!) and significantly alter the naivě prcdic- 
tions for the radiative widths. 70 Quantitative esti- 
mates are needed. 

It will also be useful to resolve the experimental 
controversy about S(^(3770) — > tt+tt- J/tp), since 
the width for i/j(1 3 D 2 ) — ► tt+tt- J/tp is related to 
the width for ^(l 3 .Di) tt+tt- J/ý. The CLEO-c 
upper limit on the ^>(3770) — ► tt+tt- J/tp branch- 
ing ratio (see the previous section) suggests that 
the latter is not necessarily much larger than the 
value induced by the 2 3 S\ — l 3 _Di mixing, while 
the BES measurement indicates a rather large rate 
for direct tttt transitions between the charmonium 
triplet ID and IS* states. The BES result implies 56 
B{ip{l 3 D 2 ) -» tt+tt- J/tP) ~ (20-40)% which would 



make it easier to accommodate BELLE's results in 
the 1 3 L> 2 intcrpretation of the X(3872). The ob- 
served dipion mass distribution does not fit the shape 
for 1D 2 — > IS transitions predicted by the multipolc 
expansion model (see Fig. 15). However, some dy- 
namical effects beyond this model can alter the dip- 
ion distribution. Finally, the íp(l D 2 ) interpretation 
is disfavored by the mass predictions from potential 
model. This is illustrated in Fig. 16. All potential 
models but one predict the ip(\ 3 D 2 ) mass to be about 
70 MeV lower then the A (3872) mass. The model 
by Fulcher 71 predicts this mass right at the observed 
value, however it overestimates the ip(l 3 Di) mass by 
similar amount. In other words, none of the exist- 
ing calculations can accommodate the X(3872) and 
-0(3770) as spin 2 and 1 members of the 1 3 -Dj triplet. 
Coupled channel effects and l 3 Z?i — 2 3 5i mixing can 
increased the mass splitting relatively to the naivě 
potential model calculations. Quantitative estimates 
of these effects are needed. 

The mass of the Xc(2 3 -Pi++) statě is predicted by 
the potential models to be higher than the X(3872) 
mass (see Fig. 16). Thus, significant coupled chan- 
nel effects would need to be at work for this inter- 
pretation as well. Decays of B+ to K+xA^- 3 'Pi) 
are observed with a rate comparable to K+tp(2S) 
and K+J/tPilS) 20 Therefore, decays to K+xd^Pi) 
should also occur. If the A (3872) is the Xc(2 3 -Pi) 
statc then the ratio £(X(3872) -> 7J/0(l 3 Si))/ 
B{X{3872) -> tt+tt- J/V>(l 3 Si)) should be large. 
The photon transition here is of an El type and is ob- 
served with a 21% branching ratio in the bb systém. 
Barnes and Godfrey predict 3.5% for the Xc(2 3 í ) i) 
statě. 56 The Xc(2 3 Pi) -> tt+tt- J/ip(l 3 Si) transition 
is of an isospin violating type and the branching ra- 
tio must be small even in the charmonium systém. 
The BELLE experiment does not see evidence for 
significant A(3872) — > jJ/tp(í 3 Si) branching ratio. 
Quantitative results should be soon available. 73 

Finally, A (3872) could be a cčg hybrid statě, 64 ' 56 
though the masses of the hybrid states are predicted 
to be significantly higher than the observed mass. 

The molecular, conventional charmonium or hy- 
brid charmonium interpretations of the X(3872) 
should not be viewed as mutually cxclusivc options. 
For examplc, the production rate for molecular char- 
monium in B decays is expected to be small. How- 
ever, the production rate can be enhanced by mixing 
of the molecular systém with conventional charmo- 
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nium states 74 ' 64 ' 65 (e.g. with Xc(2 3 í , i) for the 1++ 
molecule). Mixing would influence the pattern of de- 
cay branching ratios as well. 74 ' 64 ' 56 

More experimental studies of A (3872) produc- 
tion and decays are needed to clarify nature of this 
statě. 
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Figuře 16. The predicted masses (vertical axis) of the i/>(l 3 Di) 
(squares), ip(l^D2) (triangles) and Xc(2 3 Pi) (circles) states 
versus overall quality of the potential model expressed as RMS 
of the differences for the center-of-gravity masses of long-lived 
cč states (horizontál axis). The measured masses of i/>(3770) 
and Jí (3872) are indicated by horizontál bars. Prcdictions for 
the 1S,2S and IP masses are included in the RMS calculation 
for all models displayed here. 72 

9 Summary and Outlook 

Heavy quarkonium physics has been recently exper- 
imentally revitalized. Largo data samplcs colleeted 
in e+e~ and pp annihilation by BES-II (cč), CLEO- 
III (66) and E835 (cč) are still being analyzed. The 
_B-meson gateway to charmonium is now wide open 
with 3 x 10 8 B decays recorded by BELLE and 
BABAR. Thcre has been similar progress in theory 
by the development of NRQCD and much improved 
calculations with QCD on the lattice. 

More results for charmonium are expected in the 
future. In addition to more results from the BES- 
II experiment, the CLEO-c program will contributc 
greatly. The hrst CLEO-c results have been prc- 
sented. Run-II data from CDF are proving to be 
uscful as well. In the farther future, BTeV and LHCb 
are likely to contribute to quarkonium physics. The 
BES-III/BEPC-II project was recently approved in 
China. There is also a proposal for a new dedicated 
pp machine to explore charmonium physics (PANDA 
at GSI). 

Prospects for more 66 data are not as well de- 



fined. In principle, CLEO at CESR could go back 
to the high cnergy running and aceumulate more 
statistics for the narrow T resonances. BELLE and 
BABAR could aceumulate T data with even higher 
rates if their S-factories are ever utilized to produce 
these states. 

The discovery of the A (3872) particle by 
BELLE, and its likely compound nature, underlines 
the importance of the heavy quarkonium systems. 
Previously charmonium played a crucial role in the 
solidification of the naivě quark model. It will not 
be a big surprise if it provides the first convincing 
proof for the existence of hadronic systems beyond 
the naivě quark model. 
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